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Abstract. Recurrence and explicit formulae for contractions (partial traces) 
00 , of antisymmetric and symmetric products of identical trace class operators 

C ~ ) , are derived. Contractions of product density operators of systems of identical 

C ~ ~) • fermions and bosons are proved to be asymptotically equivalent to, respec- 

ts) ' tively, antisymmetric and symmetric products of density operators of a single 

I > i particle, multiplied by a normalization integer. The asymptotic equivalence 

^ , relation is defined in terms of the thermodynamic limit of expectation val- 

1 ues of observables in the states represented by given density operators. For 

some weaker relation of asymptotic equivalence, concerning the thermody- 
l/"") namic limit of expectation values of product observables, normalized antisym- 

Os) ■ metric and symmetric products of density operators of a single particle are 

shown to be equivalent to tensor products of density operators of a single 
particle. 

This paper presents the results of a part of the author's thesis [W. Radzki, 
Kummer contractions of product density matrices of systems of n fermions 
-t— J ■ and n bosons (Polish), MS thesis, Institute of Physics, Nicolaus Copernicus 

£h ' University, Toruh, 1999]. 

=5 

£' 

1. Introduction 

CN ■ 

k>" ' This paper, presenting the results of a part of the author's thesis [10], deals 

with contractions (partial traces) of antisymmetric and symmetric product density 
operators representing mixed states of systems of identical noninteracting fermions 
■ and bosons, respectively. 

| If H is a separable Hilbert space of a single fermion (boson) then the space of the 

. n-fermion (resp. n-boson) system is the antisymmetric (resp. symmetric) subspace 

' H An (resp. Ti. Vn ) of Ti® n . Density operators of n- fermion (resp. n-boson) systems 

. are identified with those defined on H®" and concentrated on 7i A " (resp. TL Vn ). 

Recall that the expectation value of an observable represented by a bounded 
selfadjoint operator B on given Hilbert space in a state described by a density 
operator p equals TiBp. If B is an unbounded selfadjoint operator on a dense 
subspace of given Hilbert space, instead of B one can consider its spectral measure 
£b(A) (which is a bounded operator) of a Borel subset A of the spectrum of 
B. Then Tr Eb(A)p is the probability that the result of the measurement of the 
observable in question belongs to A [9]. 

fc-particle observables of n-fermion and n-boson systems (k < n) are represented, 
respectively, by operators of the form 



X 



r£ b = a 



{ £> (B ® A%\ Ft B = S% } (B <g> S% ] (1) 

(multiplied by (")), where and S ( ^ ] are projectors oiU® n onto Ji Kn and H Vn , 
respectively, / is the identity operator on TL and B is a selfadjoint operator on 
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W® fe (see [4]). Operators (1) are called antisymmetric and symmetric expansions of 
B. In view of the earlier remark it is assumed that B is bounded. The expecta- 

A V 

tion values of observables represented by rj? B and T£ B in states represented by 
n-fermion and n-boson density operators K and G, respectively, can be expressed 



Tr K r]J B = Tr BL*K, Tr G T% B = Tr BL^G (2) 

(see [4, Eqs. (1.7), (3.19)]), where fc-particle density operators L^K and L^G are 
(n,k)- contractions of K and G (see Definition 2.1), called also reduced density 
operators. Such operators were investigated by Coleman [1], Garrod and Percus [2], 
and Kummer [4]. 

In the present paper particular interest is taken in the case when K and G are 
product density operators, i.e. they are of the form 

K = -J—/", G = —^—p Vn , (3) 
Tr/9 A » P Trp v " F ' V ; 

where p An = A% ) p® n A% ) , p v " = sff p^sff , and p is a density operator of 
a single fermion or boson, respectively. The first objective of this paper is to find the 
recurrence and explicit formulae for L n K and L^G for K and G being, respectively, 
antisymmetric and symmetric products of identical trace class operators, including 
operators (3). The explicit form of the operators L^K and L^G proves to be quite 
complex. However, they can be replaced by operators with simpler structure if only 
the limiting values of expectations (2), in the sense of the thermodynamic limit, 
are of interest. The second objective of this paper is to find that simpler forms 
of contractions L^K and L^G for product density operators (3), equivalent to the 
complete expressions in the thermodynamic limit. 

The problems described above have been solved for k = 1, 2 by Kossakowski and 
Mackowiak [3], and Mackowiak [(>]. The formulae they derived were exploited in 
calculations of the free energy density of large interacting n-fermion and n-boson 
systems [3, 6], as well as in the perturbation expansion of the free energy density 
for the M-impurity Kondo Hamiltonian [8]. In the case of investigation of many- 
particle interactions of higher order [15, 14, 7, 13], or higher order perturbation 
expansion terms of the free energy density, the expressions for (Tr p An )^ 1 L^ l p An 
and (Tr / 9 Vn )~ 1 L^p Vn with k > 3 prove to be needed in the canonical and grand 
canonical ensemble approach, which is the physical motivation for the present paper. 

The main results of this paper are Theorems 3.1, 3.4, 4.9, and 4.14. 



2. Preliminaries 
In this section notation and terminology are set up. 

2.1. Basic notation. Let (TC, (•, •)) be a separable Hilbert space over C or R. The 
following notation is used in the sequel. 
/ - the identity operator on Ti., 

B(Ti.) - the space of bounded linear operators on 7i with the operator norm ||-||, 
T(TL) - the space of trace class operators on H with the trace norm Tr |-| , 
B*(TL) - the space of bounded sclfadjoint operators on H, 
Bt >0 (H) - the set of nonnegative definite bounded selfadjoint operators on Tt, 
V(TL) - the set of density operators (matrices) on H, i.e. 

V(H) = {D e T(7i) | D = D*,D>0,TrD = l}. 
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Set H® n = T-C ig) ■ • • ig) "H and denote by S n the group of permutations of the set 

n 

{1, . . . , n} . Let A%\ S$ £ B{H® n ) be the projectors such that 
A$(ipi <8> •• • ® ip n ) = —. Sgn7r^(i) <8> • • • ® ^(n), 

for every ^i, • • • , ip n G The closed linear subspaces H An = A { ™ ] H® n and H yn = 
S^TL® n ofH.® n are called the antisymmetric and symmetric subspace, respectively. 
The antisymmetric and symmetric product of operators B E B (Tt® k ) , C £ 

£(W® m ) are defined as B A C = A ( £ +m) (B g> C)A^ +m) and BVC = S% +m) (B ® 
^^(fe+m)^ respectively _ It is assume d #An = BA. . .AB y , B yn = B V . . . V £ , 

n n 

and B A1 = B vl = £?. Clearly, if B £ then B An = A { £> B® n = B®"^ 3 , 

B v " = S$B® n = B® n S%\ and if B £ S^ (W) then B An ,B Vn £ ^(H®"). 

Set R + = [0, +oo) and R + = R + U {+00} . The product of measures fi, /ii is 
denoted by /i ® /^i and /i®" stands for /i (g> • • • <E> [i . In subsequent sections use is 

n 

made of product integral kernels, described in Appendix A. 



2.2. Contractions of operators. The definition and basic properties of contrac- 
tions of operators are now recalled for the reader's convenience. A discussion of 
these properties was carried out by Kummer [4, 5]. 

Let Tt be a separable Hilbert space over the field K = C or R. 

Definition 2.1. Let k, n £ N, k < n, and K £ T(H® n ). Then the {n, k)- contraction 
of K is the operator L*K £ T(TL®) such that 

VceB(W®*): Tr H ®4C®I®( n - k) )K = Tr n ® k CL*K. (4) 
It is also assumed L"_ftT = K. 

Remark 2.2. The operator L^K always exists and is defined uniquely by Eq. (4). 
Ij^K is a partial trace of K with respect to the component 7-^®(™-' c ) of 7i®" = 
n ®k (g n ®(n-k)^ \iu = H Y ■= L 2 (Y,n), where the measure /x is separable and 
(7-finite, and /C is a product integral kernel of K (see Appendix A) then L^K 
has an integral kernel /Co given by formula (60), according to Lemma A. 5 and 
Corollary A. 6. 

Under the assumptions of Definition 2.1 one has Tr^kL^K — Tr H (s^K, and 
if p £ N, k < p < n, then h k (L%K) = h k K. Moreover, if if € B*(H® n ) then 
L k K £ B*{H® k ), and if K £ B*> {H® n ) then L*K £ B*> (H® k ). 

Contractions of density operators are called reduced density operators. Contrac- 
tions preserve the Fermi and the Bose-Einstein statistics of the contracted oper- 
ator, i.e. for K £ A^ ) T(n® n )A { ^ ) and G £ s£ } T(H® n )S% one has L k K £ 
A$T{-H® k )A$ and L k G £ S^ ) T(n® k )S i ^ ) . For such /f and G Eqs. (2) hold. 

The following theorem is a part of Coleman's theorem [1, 4]. 

Theorem 2.3. Let n £ N, n > 2. For every (n-fermion) density operator D £ 
V (H® n ) , D = A { £ } DA { £\ one has \\L^D\\ < ± ||D|| . 
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3. Recurrence and explicit formulae for contractions of products 

of trace class operators 

In this section recurrence and explicit formulae for contractions of antisymmetric 
and symmetric powers of single particle operators are derived. 

In the whole section use is made of the Hilbcrt space Hy '■= L 2 (Y,fi) over the 
field K = C or R, where the measure p is separable and er-finitc. 

The following theorem in the case of k = 1,2 was proved in [3, 6]. 

Theorem 3.1 (Recurrence formulae). Let p G T(Hy)- If k,n S N, 1 < k < n, 

then 



©iV" = Gift (t!/ K1) )Ap 

- (V) (Lt lP « n -») (i^-V ® p) A™ , (5) 

(fe) L nP V " = (fell) (itV^Vp 

+ (l*.^"- 1 )) (Z®^- 1 ) ® p) S« , (6) 
nL„p A " = (Trp^"- 1 )) p- (n - 1) (l^p^"" 1 )) p, (7) 



nL„p v " = (Trp^- 1 )) p+ („ - 1) (L^p^- 1 )) p. (8) 

Proof. Let 7 x V ^ I be a product integral kernel of p. For every m G N define 
the mapping g Am : Y m x Y m — > K by the formula 

■■• g{xi,y m ) 

Q{x m ,Vi) ■ ■ ■ g(x m ,y m ) 
Then the mapping K, : Y n x Y"™ — > K given by 

1^/ \ 1 An i ^1 5 ■ ■ ■ 7 \ 

n\ \ yx, ■ ■ ■ ,y n J 

is a product integral kernel of p A ™ = A^' p® n . 

Eq. (5) will be first proved for n > fc + 1. In view of Remark 2.2, an integral 
kernel C : Y k x Y k -> K of (£)L„p An can be given by 



Am ( ] =det 

V yi,---,y m J 



C(x',y') = ( n k ) / )C(x',x",y',x")d^ n -Q(x") 



for /z® 2fe -a.a. (a;', y') G Y" fc x Y~ fc . Performing fc! permutations of the first fc rows and 
A;! permutations of the first k columns of the determinant defining /C and expanding 
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that determinant with respect to the fcth column one obtains 
C(xi, ...,x k ,yx,.. . ,y k ) 

= (1)37^2 J2 sgnTrsgncr^-l)^' f Q{x 1x(j)l y <k) ) 

■ n A (™ — ( Xw (^) ' ' ' ' ' X ir(J—l) ' x ir(j+l) j ■ ■ ■ 1 x ir(k) 1 x k+l 1 ■ • ■ 1 %r. 
V Vail), ■ ■ ■ ,ya(k-l),X k+1 , . . . ,X n 

dfx^ n - k \x k+1 ,...,x n ) 

,„s 1 1 



sgnTrsgncr Y (-l) fe+J / Q(xj,Va(k)) 

- e ■ . . -1 JY ri ~ k 



ykl n\(k\) 2 

v ' ir,a£S k j=k+l 

■ g) A (" _1 ) [ Xn ( 1 ) : ' ' ' ' X n(k)i ^fc+l, ■ ■ ■ 1 x j-li Xj+1, ■ ■ ■ , X r , 
\ Ucr(l), ■ ■ ■ ,V(T(k-l)> x k+l, ■■-,X„ 

dn^ n - k \x k+1 ,...,x n ). (9) 

Consider the first term on the r.h.s. of Eq. (9). In all summands of Sj=i except 
the last one the (k — l)th row of the determinant (containing the variable aWfc)) 
can be shifted into the jth position, changing thereby the sign of the determinant 
by (— l)( fc_2 ) — Ci — — (_i)-fc-j+i. Then the first term of sum (9) assumes the form 

1 1 



J2 sgn^sgn ( 7^(-l) fe ^(-l) k 3+1 J nk Q(x A3) ,y a[k) ) 

A(n— 1) [ x tt(1) 7 ■ • ■ i X 7r(j — 1) ) ^7r(fc) ; ^7r(j+l) j * * • ) ^7r(fc — 1) j j • * • ; -^r 

\ 2/<t(1)> ■ ■ ■ iy<7(fc-l),Zfc+l, ■ ■ ■ ,Xn 

d^ n - k \x k+1 ,...,x n ) 
+ ^^!(fcl)2 51 sgn 7 rsgn ( r(-l) fc+fe ^^^(a; 7r(fe) ,j/ CT(fe) ) 

. A(n-1) j x ii(1)' ■ ■ ■ > ;r 7r(fc-l)) :c fc+l) ■ ■ ■ f^ra \ ^®(n-k) ( Xkl x n ) 
\ 2/cr(l), • • • > J/<r(fc-l) > ^fc+l > • • • > x n J 

(10) 

Let Tj k S S'fc denote the transposition j <-> fc for j < fc (then (-l) fe +.J(-l)- fe -j+i = 
( — 1) = sgnTj/j) and the identity permutation for j = k (with sgnT^^ = 1). Ex- 
pression (10) can be written as 

k 1 1 f 
HCfc)ri7Tn2 (sgnTTSgnTjfcJsgno- / g(a; (7roTjfc)(fc) , y a(fc) ) 



uv n! (fc!) 2 - 

n A(n-l) ( E(noT jh )(l)> ■ ■ ■ , x {noT jk )(k-l), x k+l, ■ ■ ■ ,X n 
Vcr(l), ■ ■ ■ ,y<j(k-l)i x k+l, ■ ■ ■ , x n 

d^ n ' k \x k+1 ,...,x n ) 



Q 



l-DjUv J2 ^rsgnag(x T{k) ,y^ k) ) j 

d^- 



Q 



A(n-l) 



x t(1): ■ 
V<r(l), ■ 



> X T(k-l)i x k+l, ■ 
,y<r(k-l), x k+l, ■ 



■ x n 



-1)! 

~ k \x k+ i, ...,x n ) 



(11) 
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The function V\ : Y k x Y k — > K, such that V\{x\, . . . , Xk,Vi, ■ ■ ■ , Vk) is p^ 2k -a..c. 
equal to expression (11), is an integral kernel of the operator 

which appears on the r.h.s. of Eq. (5). 

Consider now the second term of the sum on the r.h.s. of Eq. (9) . One can change 
the indices of the integral variables x^+i, ■ ■ ■ ,Xj in all summands of X}j=fc+i except 
the first one, according to the rule Xj — > Xk+i —* Xk+2 — * • • ■ - * Xj for the jth 
summand, and simultaneously change the order of the columns of the determinant 
inversely (which changes the sign by (— l)^ -1 ) - * = ( — l)( fe+1 )~ : '). The resulting 
sum X]j=fc+i then contains n — k terms identical to the one with j = k + 1. Thus 
the second term of sum (9) equals 

- ( n - fc ) (fe) 3 7Zu2 E Sgn 71 Sgn ° i g( - Xk+1 ' ) 

1 j 7r,(7 6 S fc JY "~ k 
A(n-l) / ^(l)! ■ • ■ j ^^(fc) ? :E fc+2, ■ ■ • ,X n \ i ®(n-fe) 



2/(7(1) J • ■ ■ >V(r(k-l)> X k+l, ■ ■ - ,x n 

• p A< ' n_1 '' [ 2-1; ■ ■ ■ ; x fc; a-fc+Z; ■ • • j 

V 2/(t(1); ■ ■ ■ ! 2/<T(fc-l)> £fc+l> ■ ■ ■ i^n 
d/i ®(n-i-fc) (a , fc+2; _ _ : ^ d/i (a; fc+1 ). (12) 

The function P 2 : ^ fc x ^ fc — * K, such that Vi{x\, . . . , Xk, Ui, ■ ■ ■ , Uk) is /i® 2fe -a.e. 
equal to expression (12), is an integral kernel of the operator 



-(V) (lJU^™- 1 )) (l® {k -^®p) A 



(fe) 



which occurs on the r.h.s. of Eq. (5). One concludes that the kernel C of the 
operator on the l.h.s. of Eq. (5) is /x® 2fe -a.e. equal to the kernel V\ + V2 of the 
operator on the r.h.s. of Eq. (5), which proves the equality of both operators. 

The proof of Eq. (5) for n = k + 1 and the proof of Eq. (7) proceed analogously. 

Similarly, the proof of Eqs. (6), (8) is accomplished by changing the product A 
into V and replacing determinants in all formulae by pcrnamcnts, defined for every 
complex matrix A = [ajj]™ =1 as 

perA = E a 7r(l),l • ■ ' 0>ir{m),m- 

Notice that signs of permutations are omitted in this case, similarly as the multi- 
pliers ±1 in the Laplace expansions. □ 

Lemma 3.2. Let fc, to e N, 1 < k < m, p <ET~ (Hy) , jk £ {k, . . . , m} , and 



ifc-i jk-1-1 32-1 

E E 

jk-i=k — l jk-2 = k-2 31=1 



R:= ) > ...) fP 1 ®p> 2 - jl p">-^-i 



(for k = 2 the only summation index is jk—1 = ji)- Then A^ y R = RA^ and 
R = RS& ■ 

The proof of the above lemma consists in demonstrating the invariance of R 
under permutations of factors in the tensor products. To this end it suffices to 
observe that R is invariant under transpositions of neighbouring factors. 
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Lemma 3.3. Let p G T(Hy), £ A ■= Trp Aa , g := Tr p Vs for s E N, Co : = ^ 
Cq := 1, and 

m jp-l 32 - 1 

K p (p) ■■= E E • • • E ^- jp (- 1 ) p+jp p jl A ^ 2 " jl A ■ ■ • A p**"*- 1 . 

3p=p jp-i=p— i ii =1 

m i P - 1 J2 - 1 

re (p) ; = E E • ■ ■ E v ^ a_il v • • • v p 1 ''*'- 1 

3p=p jr-i=p—i ii =1 

for p, m € N, p < m. (For p = 1 t/ie only summation index is j\ and the summation 
runs over the operators p 31 .) If '2 < p < m then 

nt p (p) = (K^ip)) a P {KIM) (i mp - 1] ® P )A ( t (13) 

and 

< p (p) = (K^ip^p+iKiM) (i^- 1) ^p)s^ y . (w) 

Proof. Eq. (13) will be first proved for p > 2. One has 

m — 1 tja— 1 _ 1 £2—1 

n^( P ) = C-pPA...Ap+E E E - E 

lp=P lp-l=p — 1 ip-2=P — 2 (l = l 

£ A l _i I ,_l(-l) p+i, ' + V 1 A p' 2 "' 1 A ... A ph-i-h-2 A p'p-'p-i+i 

m — 1 Jp — 1 ^p — 1 — 1 ? 2 — 1 

= C- P pA...Ap+E E E • E 

(p=P Jp-i=p— 1 ip-2=p— 2 (i = l 

C-/ p -i(-l) P+ ' P+ V Zl A p' 2 ~ Zl A ... A p^-i-'p^ A p'f-'f"^ 1 
m— 1 ip — 1 ip-2— 1 Z2 — 1 

+ E E E -E 

ip=P ip-2=p-2 ip_3=p — 3 il=l 

C-/p-i(-l) P+ ' P+ V Zl A p' 2 ~ Zl A ... A p'p-^-'p-a A p h~h-2 A p ( 15 ) 
The first and the third term after the last of equalities (15) yield 
m— 1 jp-i— i 32—1 

E E -E 

jp-l=P-l jp-2=p-2 J 1= l 

A p^' 1 A ... A pfc-i-fc-») A p (16) 

for Z p = jp_i, Zp-2 = ip-2, ■ • • , h = ji- By Lemma 3.2, the second term after the 
last of equalities (15) equals 

The sum of expressions (16) and (17) is equal to the r.h.s. of Eq. (13) for p > 2. 
After simplifications the proof also applies to the case of p = 2. 

The proof of Eq. (14) is analogous to that of Eq. (13). □ 

The next theorem provides the explicit form of (n, fc)-contractions of product 
operators. The proof for k = 1, 2 was given in [3, 6]. The author of [6] emphasized 
that formula (18) for k = 2 was derived by S. Pruski in 1978. 
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Theorem 3.4 (Explicit formulae). Let k,n G N, k < n, p G T (H Y ) , £ A : = Trp As , 
4 V := Trp Vs for s G N, and £ A := 1, $ := 1. T/ien 

n 32-1 

= E E ■ ■ ■ E (-i^+V 1 a p ,2 - ii A ... A p"- 

jk j k _ 1= fc-l Jl = l 

n — (fc— 1) n— ii — (fc— 2) n—i\ 2^_ 2 — 1 n — ii ifc— l 

- E E ■■■ E E 

n = l «2=1 ife_i=l ifc=l 

e !1 -..- !t (-l) Hil+ - +1 V I1 A...A^ (18) 

and 

(n\r k Vn 

n jfc-1 J2 — 1 

= E E ■ ■ ■ E C-j-y 1 v / ,/: ' 1 v • • • v p"- 1 

n— (fc— 1) n— ii— (fc — 2) n — i"i 2 fc _ 2 — 1 n — ifc— l 

= E E ■■■ E E 

ii=l '2=1 «k-i=l «fc=l 

^-...-^V-.-Vp 1 '*. (19) 

(7*br fc = 1 £/ie onZ?/ summation indices are j\ and i\ and the summation runs over 
the operators p 31 and p n , respectively.) 

Proof. For every p, m G N, p < m, let IT A f(p) be defined as in Lemma 3.3. Then 
the first of equalities (18) can be written as 

(l)L k n p An = K k (p). (20) 

The proof of Eq. (20) will be carried out by (double) induction with respect to k 
and, for fixed fc, with respect to n > fc. 

1°. (fc = 1) This part of the proof is by induction with respect to n > 1. 

a) (n = 2) According to Theorem 3.1, 2L^p A2 = (Trp) p - p 1 = n A1 (p). 

b) Assuming validity of formula (20) (with fc = 1) for n G {2, . . . ,m — 1} , where 
m G N, m > 2, its validity will be shown for n = m. 

One has 

rn 

Kn(p) = C-ip+ E = C-ip- (n^_x(p)) P . 

Thus, according to the inductive hypothesis for n G {2, . . . , m — 1} , 

n^(p) = - (m - i) (hl^p^- 1 )) p, 

which, in view of Theorem 3.1, yields (7)L^p Am = n A1 (p). 

2°. Assuming validity of formula (20) for fc G {1, . . . ,p — 1} (and every n > fc), 
where p G N, p > 1, its validity will be shown for fc = p. For arbitrarily fixed p the 
proof will be carried out by induction with respect to n > p. 

a) (n = p + 1) By the inductive hypothesis with respect to fc and Lemma 3.3, 

KIM - (K-i 1 ) te-iP A((p+1H) ) ap-p^i^- 1 ) ®p)4£, 

hence ( p + 1 )L^ +1 p A ( p+1) = n A ^(p), according to Theorem 3.1. 

b) Assuming validity of formula (20) for n G {p + 1, . . . , m — 1} , where fc = p, 
m G N, m > p + 1, its validity will be shown for n = m. 
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By the inductive hypothesis for k € {1, . . . ,p — 1} and Lemma 3.3 one has 
According to the inductive hypothesis for n & {p + I, . . . ,m — 1} one thus obtains 

- (V) (L p m -iP A(m ^ 1) ) (^ {p - 1) ® p) 4t 

which, in view of Theorem 3.1, yields (™)L^p Aro = H$(p). This completes the 
inductive proof for Eq. (20) with respect to n > p and with respect to k. 

Now turn to the second of equalities (18). For k = 1 it is identity. Let k > 2. 

Setting ji = ii, j 2 = h + ia, • • • , jk = h + h ik or, equivalently, h = ji, 

i 2 = ]2 - ji, h = 33 - h, ik = jk - jk-i, one checks that both sides of the 
equality in question are equal to 

n-(fe-l) n-(k-2) n -l n 

E E ■■■ E E 

ii = l 3a=ji + l jfc-i=jfc-2 + l Jfc=Jfe-i + l 

£ A _ ifc (-l) fc+ifc p Jl A p n ~ jl A ... A p»'*-J*-i . 
The proof of Eq. (19) is analogous to that of Eq. (18). □ 



4. Asymptotic form for contractions of product states 

The explicit forms of the contractions of product states given by Theorem 3.4 
are quite complex. In the present section they are replaced by simpler operators, 
equivalent in the thermodynamic limit. The main results in this section are Theo- 
rems 4.9 and 4.14. 

In what follows use is made of the Hilbert space Tin '■= L 2 (Vl,p) (over C or 1), 
where the measure p is separable, er-finite, and satisfies the condition /i(f2) = +oo. 
For every //-measurable subset Y C f2 it is assumed Hy L 2 (Y,p). 

Let AA(£l) be a fixed family of measurable subsets of f2 such that < fJ,(Y) < +oo 
for every Y G (it can be the family of all such subsets). Fix d € K, d > 0, 

and assume that there exists a sequence {y n }„ 6N C A4(Cl) such that >y , — > d as 
n — > oo. 

Definition 4.1. Fix d G M, d > 0, and let {b Ytn }, Y n)eM(Q)xN ^ e a f am ily of 
complex numbers. A complex number b is said to be the thermodynamic limit of 

Ti 

this family if for every sequence {Y n \ (=M C M.(d) such that lim — — - = d the 

n->oo fX(Y n ) 

condition lim by n n = b is fulfilled. In such a case 6 is denoted by d — lim 6y n . 

Special attention will be given to the families of complex numbers of the form 
Tr {h k n K Yn )C Y , where k, n G N, n > k, K Y , n G T{Hf l ), and C y € k ). 

Definition 4.1 does not guarantee the convergence of families {&y>} of interest 
in physics. To obtain such a convergence, additional conditions (such as conditions 
of uniform growth [11]) are usually imposed on the sequence {^} ngN in question. 
However, those additional conditions do not affect considerations in this paper. 

Expression of expectation values of observables in mixed states by using trace, 
mentioned in Introduction, is the motivation for the following definition. 

Definition 4.2. Fix k G N and d G R, d > 0. Families {^- Y n} ^ Y n )£M(n)xN an ^ 
{-By,n}(y n )eM(Q)xN °^ °P era tors A Yni B Y n G T(Tiy k ) are said to be asymptotically 
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equivalent (symbolically: A Y<n ~ B Yn ), if for every family {C-Y,n}^Yn)eM((i)xN 

of 

operators C Yn S B{TLy k ) with uniformly bounded operator norms one has 

d - lim Tr (A Y , n - B Y , n )C Yn = 0. (21) 

n,/^(y) — >oo 

Condition (21) is required to hold in particular for families {(^Y,n}(Yn)eM(Q)xti 
such that Cy. n = Cy,m for all Y E M{Vl), n, m E N. 

Remark 4.3. The authors of [3, (>] used some different definition of asymptotic 
equivalence of families of operators, closer to Definition 4.10 in this paper. 

Remark 4.4. For fixed k E N and d E R, d > 0, the relation « is an equivalence 
relation. If A Y ^ n ~ B Yn then for every family of operators Cy^ n as in Definition 4.2 
the limit d — lim Tr Ay in C Yn exists iff the limit d — lim Tr B Y ^ n C Yn exists, in 

n,/^(y) — >oo n,/i(y) — >oo 

which case both limits are equal. Notice also that if A Y ^ n s» B Y ^ n then A Yn +D Yn 
By.n + D Yn and aAy.„ ~ aB Yn for every family {-Dy^nl^^gju^xN 
and a G C. Furthermore, for every family {Ay in } ^ n )g,vi(Q)xN 

C T (W® fe ) with 

uniformly bounded trace norms Tr |^4y n | and for every sequence {a n } ngN C C 
convergent to a E C one has a n A YjTl w aA Yn . 

Lemma 4.5. Let {^y,n} (y ,„ )6jV((n)xN and {5y,4(y,„) e M(n)xN 6e as in De/?m- 
iion ^.,2. Then 

d - lim |Ay,n - By,n\ = Ay,„ By, n . (22) 

n,^i(y) — >oc 

Moreover, if the operators A Yn , B Yn are selfadjoint then 

A Yn ^B Yn d-lim \A Yn - B Yn \ = 0. (23) 

n,/^(y) — >oo 

Proof. Implication (22) follows from Definition 4.2 and the estimate 

|Tr {Ay, n - By n )Cy, n \ < \\Cy, n \\ Tr \Ay n - B Y , n \ . 

Now assume that A Y>n m By n , which is equivalent to the condition 

D Y>n » 0, (24) 
where Z?y,n := Ay„ — By^. The operators -Dy,n have the spectral representations 



oo 



£>y n = Ai(F, n)P ipi( y„ ) , 
1=1 

where P lfii (Y,n) are the projectors onto orthogonal one dimensional subspaces of 
eigenvectors ipi(Y,n) of D Y ^ n , corresponding to eigenvalues Xi(Y,n) E R. Since 

oo 

\\(Y,n)\ = Tr \Dy n \ < +oo, for every (Y, n) E A4(Cl) xN there exists m(Y", n) E 

i=l 

oo ^ 

N such that } |Aj(y,7i)| < — . Thus the operators 

n 

i-m(y,n)+l 
m(y n) 

fy,n= X! M Y > n ) P <Pi(.Y,n) 
i=l 

satisfy the condition 

oo 

rf-lim Tr \D Y n -F Yn \ = d - lim ^ \Xi(Y,n)\ = 0, (25) 

„, M (y)^oc n ^( y ^°°i= m (y,r l )+i 
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which, in view of implication (22) proved and condition (24), yields F Yn ps D Yn ps 
0. In particular, 

d - lim Tr Fy n Cy„ = 0, (26) 

Tl,Lt(Y )— »oc 

where 

m(y,n) 

Cy,n = X! ^(^(^i 71 ))^^.")' \\ C Y,n\\ = 1- 

i=l 

Observe that Tr i*y >n Cy )n = Tr |Fy,„| , hence condition (26) gives 

d-lim Tr \F Y , n \ = 0. (27) 

n,/j,(y)— >oo 

Since Tr |Dy, n | < Tr \D Y . n - F Yn \ + Tr \F Y ,„\ , conditions (25) and (27) yield 
d — lim Tr \A Y ^ n — By, n \ = d — lim Tr |-Dy n | = 0, 

n,/^(y) — >oo n.^{Y) — >oc 

which proves implication (23). □ 

The following lemma follows from Lemma 4.5. 

Lemma 4.6. Fix k,m G N. Let {A Y , n } {Ytn]eM{n)xN and { B Y,n} {Y ^ )£M(n)xN be 
families of self adjoint operators A Yn ,B Yn G T(7iy fe ) such that A Yn ps B Yn , and 
let {D Y ,n}( Y ,n)eM(Q)xN be a f am ^V °f operators D Y>n G T(TL Y m ) with uniformly 
bounded trace norms Tr |-Dy n | . Then 

Ay tn <g> L>y,n W Sy n (g) Dy >n , Dy,„ (g) Ay,„ PS L>y,n ® By,n, 
A Y: n A £>y„ « Ry,n A Dy )?l , £>y n A Ay, n ps Dy, n A -By, n , 
Ay,„ V £>y n « Ry,n V D Y>n , D Y ^ n V j4y,„ ps Dy,„ V Ry, n - 

In the sequel {py}y e x(j2) denotes a family of nonnegative definite selfadjoint 
operators p Y G T (Wy) , and for every (Y, n) G .M(f2) X N it is assumed that 

tA 1 tV 1 Al VI 

?y,o •= J-i ?y,o ■= x > •= P^j Py ■= Py, 



& >n :=Trpp>0, &, n :=Trp?>0, 

fV 

»y,n-l 



sy.ra-l v ?1 



a y,n — *A ' °Y,n '~ tV 

sy,n ?y,n 
The objective of this section is to find density operators of the most simple form 
which are asymptotically equivalent to the operators 

A W _ T fe ( 1 An \ V(*) _ Tk ( I Vn \ 
a Y,n— L n I ^\ Py I I °Y,ri — L n I Py I ' 

defined for fixed fc G N and every (F, n) G .M(fi) x N, n > k. 

Remark 4.7. For every (F, ?i) G »A//(£1) x M the operator / -4- s Y n j r \PY is invertiblc 
and || (/ + s Yn+1 p Y y 1 \\ = 1. Furthermore, if s y „ +1 ||py|| < 1 then I - s Yn+1 p Y 
is invertiblc and ||(J - Sy )n+1 py) _1 || = (1 - Sy„+i Hpyll)" 1 - 

The next theorem is a version of a theorem studied in [3, 6] (see Remark 4.3). 

A (i) A (i) 

Theorem 4.8. If 'ffy„Wy n+ i and ffte rea/s Sy ra+1 ||py|| , (Y, n) G .M(fi) x N, are 
uniformly bounded then 
A (i) 

°Y,n (I + 4,n+iPY) « (n + l) -1 < B+ iPy, (28) 
£jLra (n + l)- 1 4„ +1 py(/ + Sy.n+iPr)" 1 . (29) 
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v (l) v (l) 

V a Yn~ cr Yn+i an d ^ e r eals s Yn \\py\\ , (Y,n) € A4(Q) x N, are uniformly bounded 
then 



/(i) 



a Y.n {I ~ Sy.n+lPY) « (ft + 1) S Yn+lPY - 



(30) 



J/, additionally, s Yn \\py\\ < e /or some e < 1 and every (Y, n) G AA{£l) x N /j/ien 



-1 „V 



iPy{I - Syu+iPy) 1 - 



A (l) At 1 ) 

Proof. By Theorem 3.1 and the assumption a Y n ~<J Y one has 



Since Tr 



v(l) 



(1) 



(n + 1) ^^n+ipy « -(n + 1) X n cr y „ (sy„ +1 py) 



(31) 



(32) 



0" 



i'.ll 



= s Yn+i \\py\\ , relation (32) 



< 4, n+ illHI Tr 

yields (28), in view of Remark 4.4. 

Now turn to the proof of relation (29). According to Remark 4.7, 



Tr 



oy -(n + 1) 1 s Y +1 py(I + s Y +1 p Y y 



Tr 



^y.n U + s y ,u+iPy) - (n + 1) l s Yn+1 p Y 



= Tr 



a Y,n ( J + Sy.n+lP^) s Y,n+lPY 



(33) 



a(!) a(!) 

The explicit form of cr Yn given by Theorem 3.4 shows that <J Yn commutes with 

a(!) 

I + $ Yn+ iP Y , and since both operators are selfadjoint, & Yn (I + s Yn+1 p Y ) is also 
selfadjoint. Thus conditions (28), (33), and Lemma 4.5 yield (29). 

The proof of relations (30), (31) runs parallel to that of (28), (29). Notice that 
in this case the expression ||(7 + s Y n+1 p Y ) _1 1| = 1 from estimate (33) is replaced 
^ IK^-Sy.n+iPy)" 1 !! = (l-Syn+ilMir 1 < (1-e)- 1 (see Remark 4.7). □ 



The following theorem for k = 2 (with the reservation of Remark 4.3) was ob- 
tained in [3, 6]. The author of [G] gave also arguments that can be used to check 
the assumptions of this theorem. 



aW a(*0 

Theorem 4.9 (Asymptotic formulae I). // fy„~ffy„ +1 for every k € N and 
s Yn \\ P y\\ < 2 for every (Y,n) € M(Q) x N 



then, for every k £ N, k > 2, 



k\ a Y n A ... A a Y ,i 

k 



v( fc ) v( fc ) 
V a Y n 7>iC7 Y n+i f or every k £ N and 

Sy n ||/9y|| < e /or some e < 1 and every (Y, n) £ A^(il) x N 

</ien ; /or every k £ N, fc > 2, 



v( fe ) 
0" 



V (l) V (!) 
y ~ fc! CT y V • • • V oy 



(34) 
(35) 

(36) 
(37) 
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Proof. First equivalence (35) will be proved. Observe that 
2Tr 



a(«) , A« a« 
°Y,n -9 ! CT y,n A ... A cr y 



Tr 



A.. .A 



< Tr 



^n) (l^- l) ®{I-S$, n+lPY ))A 
(ft -g! A ... A (^ (g-1) ® ( 7 + 4,n + iPr)) A 



(<z) 

Hy 



+ ^y.n -g! °V,n A ... A 



(?) 

Hy 



(?) 

Hy 



K" _s yn+lH Tr 



A(9) , a(!) a(1) 

°Y,n °Y,n A ... A oy 



hence 



(2-||/-<„+iPk||) Tr 



A(9) , A« A« 

°V,* -g ! ^y.n A ... A oy 



< Tr 



^ -g! C A ... A (l^" 1 ) (I + 4 n+1 py)) A« 



(38) 

Since the operators py are trace class, inf {tp,p Y ip} = 0. Thus, by assump- 

<p£H Y ; IMI=i 

tion (34) and the selfadjointncss of the operators / — s Yn+1 pY, one obtains 

II 7 - S Y.n+lPY\\ = SUp \(ip, (I - Sy +1 py)<p)\ 

IMI=i 



max { 1 - s Y n+1 inf (ip, p Y (p) , s Y +1 sup (<p,p Y (p)-l 

\\V\\=1 \\ v \\=l 

The rest of the proof of (35) is by induction with respect to k > 2. 
1°. (k — 2) By Theorem 3.1 for n > 2 one has 

( n+ l)a ( 2 ) ^n+i = ^TT CTy ' n A ^ + } S ^'"+ lPy ) 

\(2) 



1. (39) 



(n + 1) : 



i(a) ^y[n(l®^ n+1 p Y ))A^ Y . 



(40) 



Assumption (34) gives Tr 



(2) 

*y,« (^®(4n+lPy))^ 



< 



A( 2 ) 

°y,r» 



< 



, \ . A ( 2 ) A ( 2 ) 

2, hence, by Eq. (40), Remark 4.4, and the assumption o~ Y n ~o Y „ +1 , one obtains 



£y,l + ^yi (-f ® (syn+iPy)) A^ « 2-^-j- a Yn A ((n + 1) l s Y ^ +l p Y ) 



Thus, in view of equivalence (28) from Theorem 4.8 and Lemma 4.6, one has 



d2) 



v(2) 



^y,™ + CT y™ i 1 ® (*yin+iPv)) A 



(2) 
Hy 



x (l) 



n + 1 



CTyn A oy U + Sv„ + i/Jy) 



(41) 



Furthermore, assumption (34) implies that the trace norms of the operators on the 
r.h.s of (41) are uniformly bounded. Therefore, according to Remark 4.4, 



'a ( 2 ) 
a 



n A« A(l) 
Y,n - 2 a Y.n A Py, 



(42) 
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A ( 2 ) a(!) aW a( 2 ) 

The explicit form of cr Y , Vy. n A Cy n given by Theorem 3.4 implies that cr Yn 

—2 g y n A <?Yn an d (-f <8> (1 + s Yn+1 py)) commute, which proves the selfad- 
jointncss of the operator on the l.h.s of (42). Thus conditions (42), (38) for q = 2, 
(39), and Lemma 4.5 yield relation (35) for k = 2. 

2°. Assuming validity of equivalence (35) for k G {2, . . . , q — 1} , where q 6 N, q > 2, 
its validity will be proved for k — q. 
By Theorem 3.1 for n > q one has 



1 



( a ) a Y,n+l - 



1 



(n + I)? v 9 



Assumption (34) implies 



(n + 1)« 



(43) 



Tr 



a(i) 



<2, 



a(«) a(«) 

hence, in view of Eq. (43), Remark 4.4, and the assumption a Y n ~ a Y n+ii 



« (n+ 9 j )g -l (g-l) ^ A ((^+l) _1 4n+lPr) • 

Thus, by relation (28) from Theorem 4.8, Lemma 4.6, and Remark 4.4, one has 



(?) 



+ £! n (l*<*- 1) ®{8$ <n+1 py))A 



(9) 
Wy 



(9-1) 



A oy (J + 4„ +1 ^) , (44) 



since the trace norms of the operators on the r.h.s. of (44) arc uniformly bounded, 

by assumption (34). Furthermore, in view of Lemma 4.6 and the inductive hypoth- 

aC?- 1 ) , , aW aC 1 ) 
esis (Ty n w (q - 1)! cr yn A . . . A oy n , condition (44) yields 



hence 



9-1 



a W 



(9) 



q! oy A ... A oy A oy (7 + s Y +lPY ) 



.(1) 



9-1 



A (9) 
fv, — 



V 



1 -g! ff n A...Aff„ (j®^ ® (J + 4,n + lPr)) ^ « 0. (45) 



a(«) A« a(!) 

From the explicit form of oy n , <Jy n A ... A cry n given by Theorem 3.4 one finds that 

°Y,n °Y,n A... A a Yn and (l®^- 1 ) ® (/ + s Yn+1 p Y )) A$ y commute, which 
proves the selfadjointness of the operator on the l.h.s of (45). Thus conditions (45), 
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(38), (39), and Lemma 4.5 yield a Y ,n~ 9- 'y n A ... A a Yn . Validity of relation (35) 



has been proved. 

Now turn to equivalence (37). Similarly to (38) one has 



(2- \l + s Y<n+1 pY ) Tr 



V (9) v( 1 ) v(!) 

°Y,n ~l\ °Y,n V ' ' ' V a Y,n 



< Tr 



/v(«) vt 1 ) vW \ / n ^ 



-lP^))^l 



Furthermore, according to assumption (36), 

2-||/+ 4, n+ iPy II > 2 - (1 + 4,n+i IMI) > 1 - e > 0. (46) 

The rest of the proof of (37) is by induction with respect to k > 2 and proceeds 
analogously to the proof of (35) with condition (39) replaced by (46) and the 
operators / =p s Yn+1 p Y replaced by / ± s Yn+1 p Y (inversion of signs). □ 

Theorem 4.9 allows to replace (n, fc)-contractions of antisymmetric and sym- 
metric product density operators by antisymmetric and symmetric products of 
1-particle contractions, respectively, if the number n of particles in the system 
is large. Further simplification, consisting in replacement of antisymmetric and 
symmetric products by tensor products, will be now proved possible. To this end 
weaker conditions on the asymptotic equivalence relation will be imposed. 

Definition 4.10. Fix k G N and d £ R, d > 0. Families {A Y ^ n }^ Y n ) £j vf(n) xN , 
{BY,n}(Y,n)eM(0)xti °^ °P erators A YjTl , B Y , n £ T (H Y ) are called weakly asymptot- 
ically equivalent (symbolically: A Yn ~ B Yn ), if d — lim Tr (A Yn — B Ytn )C Yn — 

n,[i{Y) — >oo 

k 

for every family {CV, n } (rn)eM(n)xN of operators of the form C Y>n = <^C^ n , 

»=i 

where C^ n £ S (Ky) (i £ {l,...,fc}, (Y,n) £ M(Ct) x N) arc operators with 
uniformly bounded operator norms. 

The relation ~ has the properties analogous to the properties of the relation ks 
from Remark 4.4. 

Definition 4.11. Let k £ N, k > 2. Fix tt £ S k . A set X C {1, . . . , k} is called 
a cyclic set of the permutation tt, if X = {li, . . . , Z 9 } for some l\, . . . , l q £ {1, . . . , k} , 
g £ {2, . . . , k} , such that n(l s ) = l s +i for every s £ {1, . . . , q — 1} , and 7r(Z 9 ) = l±. 
A singleton -{7} £ {1, . . . , k} such that ir(l) = I is also called a cyclic set of the 
permutation tt. 

Note that the set {1, . . . , k} from the above definition can be represented as the 
union of disjoint cyclic sets of tt. 

Lemma 4.12. Let k £ N, k > 2. If B^\ . . . , B^ eT(H y ) then 

k\ Tr (< ® • • • ® B« ) A™ = Bgair J! ^ ft B^\ (47) 

ireS fc 3=1 s=l 

where p(tt) £ {1, . . . , k) is the number of disjoint cyclic sets of tt, indexed by j, and 
qj denotes the number of elements ofthejth cyclic set ofir, which is {lj,i, ■ ■ ■ , lj, qj } , 
where 

7r (h,qj) =l oA and, for qj > 2, ir(l jtS ) = lj, s+ i, s = 1, . . . ,qj - 1. (48) 
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p(ir) p(tt) q. 

Clearly, qj = k and {lj, s } = {1, • • • , kj . 

j=l j=l s=l 

Proof. Let {v?i} i6N be an orthonormal basis of 7iy- One has 

p(ir) 

■®<)4i= e 

where 



fc!Tr 



(49) 



If qj > 2 for some j S {1, . . . ,p(7r)} then, by condition (48) and Parseval's formula. 

oo oo 

E •■• E 

it . —1 it . — 1 

oo oo oo 

= E E ■■■ E 



Performing successive summations one then obtains 

00/ / <Jj \ \ 4. 



The derivation of the above formula for qj = 1,2, after simplifications, proceeds 
analogously. This completes the proof of Eq. (47), in view of Eq. (49). □ 



Lemma 4.13. One has 



d — lim 

n,/^(y) — >oc 



a(i 

0~ 



Y. n 



(50) 



v (2) v (l) v (l) 

and if a Yn ~ 2 cr Yn V o~ Yn (see Theorem 4-9) then 



d — lim 

n,[i(Y) — >oo 



v(!) 



0. 



(51) 



Proof. To prove Eq. (50) it suffices to observe that, according to Theorem 2.3, 



1 1 



1 1 



<-7a-|I^"II<-7a- T ^" 



1 



Now Eq. (51) will be proved. Let {yi} igN be an orthonormal basis of TCy for 
fixed Y E M(Cl). Then 

'v(l) V (!) 



v(i) v (i) 
Tr 2 <r y> „ V a Y = 2 Tr ( a Y<n ® a Y<n S H 



(2) 



/ vW \ / v« 



'1,22 = 1 



Tr a. 



. v(!) V« 



(52) 
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, s v( 2 ) v(!) vM 

Taking into account Eq. (52), the relation cr Yn ~ 2 <J Yn V oy„, Definition 4.2 for 

„, v( 2 ) v(i) 

Cy,n = I , and the equality Tr 0y„= Tr oy n = 1, one obtains 



. v(!) v(!) 

d - lim Tr cr Yi „cry : „ 

71 , /J, ( Y ) — * OO 



= 0. 



(53) 



Furthermore, 



V (!) 
a Y,n V 



V (!) V(!) 



V (!) V« 



¥>, a Y,n a Y,n <P ) < Tr K.^y,,. 

for every 93 £ such that ||</?|| = 1, hence Eq. (53) yields Eq. (51). 



□ 



Notice that Eq. (50) can be also proved analogously to Eq. (51) under the addi- 

. , . A ( 2 ) „ AW aW 

tional assumption oy n « 2 oy n A oy n . 

The proof of the next theorem for k = 2 was given in [3, 6]. 



Theorem 4.14 (Asymptotic formulae II). Let k G N, k > 2. One has 

a(i) 



a(D 

k\ <J Y n A ... A <7 F 



a(i) 



a(i) 

(7 



y.ru 



(54) 



and ifd — lim 

n,/i(V) — *oo 



('see Lemma 4-13) then 



„v(i) v(l) 

ft! £Ty V • • • V a Y,n 



v(i) 



Proof. First Eq. (54) will be proved. Fix a family {Cy,„} (y ^ )£jV|(si)xN 
such as in Definition 4.10 and set 



(55) 

of operators 



B (r) A« Mr) 



r = l. 



Then, by Lemma 4.12, one has 

, /aW a(i) 



a 



CD 

Y,n 



<g>C 



0) 
Y,n 



k\ Tr LB 



(i) 



p(7r) 



TTSSfc 3 = 1 s= 



ft*,*) 



Thus 



Y,n 



aW 
a 



'Y,n 
p (tt) 9 3 - 



E 8 *™ II ^11*8 

3=1 8=1 



7reS fc 
TT^Id 



a(i) 



a(i) 

T Y. 



Ad) 



(fe) 

Y,n 



(i:) 



p(7r) qj 

E sgnTrnTrn^fe 

TTSSfc 3=1 8=1 



(56) 
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Now, let tt e Sk, TT ^ Id, be fixed. If qj = 1 for some j 6 {1,... ,p(^)} then 



Tr Y[B 



Y,n 



whereas if qj > 2 then 



y.n 



< 



< 



< 



TrB 



g 3 --l 



< 



s=l 



y,n 



a 



Tr 



Tr 



a 



(h,i) 



D 



a(i) 



a(i) 



-1 



n r 



. s=l 



y,n ' 



Tr 



a(i) 



nii^ 

s=l 



Y.n 



Since 7T 7^ Id, there exists at least one j €{!,.. . ,p(tt)} such that qj > 2, hence 



P(tt) 

J = l 8=1 



<jj — 1 \ 1i 

n n ? 

J=l 8=1 



p(tt) 

ncy n ii n 



93 -1 



and at least one exponent qj — 1 is nonzero. Thus, by the uniform boundedness 
of the norms ||Cy in || and Lemma 4.13, the termodynamic limit of the l.h.s of (56) 
equals 0, which proves the validity of relation (54). 

The proof of relation (55), after discarding the permutation signs and replacing 
A by V, proceeds analogously. □ 



Appendix A. Product integral kernels of trace class operators 

In this section theorems concerning product integral kernels, exploited in Sec- 
tion 3, are formulated. 

Fix the Hilbert space TCy L 2 (Y,fi) over the field K = C or K, where the 
measure /i is separable and cr-finite. For every weft the space 7iy™ is identified 
with L 2 (Y n , fi® n ). Unless otherwise stated, elements of L 2 spaces are identified 
with their representatives and denoted by the same symbols. 

Let K. G L 2 (Y 2 ,[i® 2 ). In the case of the integral operator K : Hy — > Hy defined 
for every <p £ Hy and ^-a.a. x € Y by 

{K<p){x) = J K{x, y)<p{y) dMv) (57) 

both K, regarded as an element of L 2 (Y 2 , n® 2 ) as well as its arbitrary represen- 
tative is called an integral kernel of K. The kernel fC is unique as an clement of 
L 2 (Y 2 ,/j,® 2 ) but a representative of K, of a special form, given in Lemma A. 3 and 
Definition A. 4, is useful in computations of the trace of K. 

Let Ti.S(TCy) be the space of Hilbert-Schmidt operators on Hy with the inner 
product defined by (A, B)-hs(Hy) := TrA*.B and the induced norm denoted by 
II'IIhschv) ■ I n t ne sequel use is made of the following theorem, the proof of which 
can be found in [12]. 

Theorem A.l. An operator K E B(Hy) is Hilbert-Schmidt iff it is an integral 
operator with an integral kernel K € L 2 (Y 2 ,[x® 2 ). Furthermore, \\K\\- hs ^ Hy ^ = 

II^-IIl^y 2 ^ 82 ) ■ 
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Corollary A.2. Let K, G £ HS(H Y ) and let K, Q £ L 2 (Y 2 , ^ 2 ) be integral kernels 
of the operators A, G, respectively. Then {K,G)-hs(Hy) = (^>^) l 2 (y 2 ^® 2 ) • 

Recall that A' £ B(H Y ) is a trace class operator iff there exist operators K\, K 2 £ 
HS(H Y ) such that K = K X K 2 . Moreover, TrA = {K{, K 2 ) hs(Hy) . This fact, 
Theorem A.l, and Corollary A. 2 imply the following lemma, in which elements of 
the L 2 space are distinguished from their representatives. The element of the L 2 
space represented by a function / is denoted by [/] . 

Lemma A.3. Let K £ T (H Y ) , K = I<iK 2 , where K 1 ,K 2 £ HS (H Y ) ■ Let 
[/Ci], [JC2] £ L 2 (Y 2 , /i® 2 ) be integral kernels of K-y, K 2 . Then for any choice of rep- 
resentatives TC\ £ [/Ci], K, 2 £ [/C2] i/ie function IC: Y x y — * K defined for ji® 2 -a.a. 
(x,y) £ y x y 6t/ 

fC(x,y) = J Ki(x, z)IC 2 (z, y) dfi(z) (58) 

is /j,® -square integrable and it is an integral kernel of A. TTie function £: Y — > K 
defined for fi-a.a. x £Y by C{x) = /C(x, a;) is fi-integrable. Moreover, 

TrA = y £(x)d£t(x) = ^ /C(x, x) d/x(x). (59) 

Definition A. 4. Under the assumptions of Lemma A. 3, the function /C given 
by formula (58) (for any choice of representatives /Ci, fC 2 of [JCi], [tC 2 ]) is called 
a product integral kernel of A. 

Notice that for /j, being the Lebesgue measure on [0, 1] x [0, 1] formula (59) is 
valid, for example, if /C is any continuous function. 

In the following lemma, which follows from Lemma A. 3, the function /Co need 
not be a product integral kernel of Kq but the integral formula for the trace of Ao 
still holds for /Co- 

Lemma A. 5. Let k, n £ N, k < n, and let IC be a product integral kernel of 
A £ T (H® n ) = T (L 2 {Y n ,n® n )) . Then the function IC :Y k x Y k -> K defined 
for n® 2k -a.a. (x',y') £ Y k x Y k by 

K (x',y')= [ JC(x\x",y',x")d^ n - k \x") (60) 

is /i® 2k -square integrable and the integral operator Kq with the kernel /Co belongs to 

iS N °f H Y 



T (Hy k ) . For every \, ip £ Hy k and every orthonormal basis {V'ijjgpj of H?^ 



one has 

00 

i=l 

The function £ : Y k — > K defined for [i® k -a.a. x' £ Y k by Cq(x') — Kq{x' ,x') is 
H® k -integrable. Moreover, 

K {x', x') d/j,® k {x') = [ Co(x') d^ k (x') = Tr A = Tr K. 
Y k J Y k 

Corollary A. 6. Under the assumptions of Lemma A. 5, if C £ B (H® k ) then 
Tr CK = Tr (C <g> I®( n ~ k ))K. 
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